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Abstract. - The tumbling dynamics of individual polymers in semidilute solution is studied by 
large-scale non-equilibrium mesoscale hydrodynamic simulations. We find that the tumbling time 
is equal to the non-equilibrium relaxation time of the polymer end-to-end distance along the flow 
direction and strongly depends on concentration. In addition, the normalized tumbling frequency 
as well as the widths of the alignment distribution functions for a given concentration dependent 
Weissenberg number exhibit a weak concentration dependence in the cross-over regime from a 
dilute to a semidilute solution. For semidilute solutions a universal behavior is obtained. This is 
a consequence of screening of hydrodynamic interactions at polymer concentrations exceeding the 
overlap concentration. 



Introduction. — Polymers in solution exposed to 
shear flow exhibit a remarkably rich dynamical behavior, 
■ as has been shown by direct observation using fluorescence 
microscopy [T]-[5]. In particular, polymers exhibit tumbling 
motion, i.e., they undergo a cyclic stretch and collapse dy- 
namics, with a characteristic frequency which depends on 
the shear rate and their internal relaxation time. This 
nonequilibrium behavior has intensively been studied for 
polymers in dilute solution [Till 3] . 

The dynamical behavior of a polymer in semidilute so- 
lution under shear flow has received far less attention 
[T4HT6] . Insight into the behavior of such systems is of 
fundamental importance in a wide spectrum of systems 
ranging from biological cells, where transport appears in 
dense environments, to turbulent drag reduction in fluid 
flow. While the dynamical behavior of polymers in di- 
lute solution is strongly affected by hydrodynamic inter- 
actions [T7HT9] . their relevance in semidilute solutions is 
less clear. 

The complex interactions in semidilute solutions ham- 
per an analytical treatment. Here, computer simulations 
are an important tool to shed light on the rich and intricate 
dynamical behavior of such systems. The large length- 
and time-scale gap between the solvent and macromolccu- 
lar degrees of freedom requires a mesoscale simulation ap- 
proach in order to assess their structural, dynamical, and 
rheological properties. We apply a hybrid simulation ap- 



proach, combining molecular dynamics simulations (MD) 
for the polymers with the multiparticlc collision dynamics 
(MPC) method describing the solvent fL8H2u] . 

By this approach, we demonstrated that polymers in 
dilute and semidilute solutions exhibit large deformations 
and a strong alignment along the flow direction in sim- 
ple shear flow |16j . More importantly, in the stationary 
state, the conformational and rheological properties for 
various concentrations are universal functions of the Weis- 
senberg number Wi c = 7t(c), where 7 is the shear rate and 
r(c) the concentration-dependent polymer end-to-end vec- 
tor relaxation time at equilibrium. Hence, with increas- 
ing concentration, hydrodynamic interactions affect the 
conformational and rheological properties only via the in- 
creasing relaxation time t(c). Experiments on DNA in 
shear flow [T5] and simulations of polymer brushes pH] 
lead to a similar conclusion. Then, the question arises 
to what extent hydrodynamic interactions are relevant in 
non-equilibrium systems. 

In this letter, results are presented for the concentration 
dependence of the non-equilibrium dynamical properties 
of polymers in shear flow by calculating tumbling times 
and oricntational distribution functions. These quanti- 
ties exhibit a dependence on hydrodynamic interactions 
in dilute solution, and are independent of such interac- 
tions in semidilute solution where hydrodynamic interac- 
tions are screened. This is supported by a comparison of 



p-1 



C.-C. Huang et al. 



non-draining and free-draining simulations. As a result, 
hydrodynamic interactions ar found to clearly contribute 
to the non-equilibrium polymer dynamics in dilute solu- 
tion beyond the change of relaxation times. 

Model and Parameters. — A solution is considered 
of N p linear flexible polymer chains embedded in an ex- 
plicit solvent. Each polymer is comprised of N m beads of 
mass M, which are connected by linear springs of equi- 
librium bond length I |16U19j . Inter- and intramolecular 
excludcd-volume interactions are taken into account by 
the repulsive, shifted and truncated Lennard- Jones poten- 
tial, with the parameter a characterizing the bead size and 
e the energy |16j . The monomer dynamics is determined 
by Newton's equations of motion, which are integrated by 
the velocity Verlet algorithm with time step h p [22] . 

The solvent is simulated by the multiparticlc collision 
dynamics (MPC) method [I5H2D] . It is composed of JV S 
point-like particles of mass m. The algorithm consists of 
alternating streaming and collision steps. In the streaming 
step, the solvent particles move ballistically for a time h. 
In the collision steps, particles are sorted into cubic cells of 
side length a and their relative velocities, with respect to 
the center-of-mass velocity of their cell, are rotated around 
a randomly oriented axis by a fixed angle a. 

The solvent-polymer coupling is achieved by taking the 
monomers into account in the collision step. To insure 
Galilean invariance, a random shift is performed at every 
collision step [23]. The collision step is a stochastic pro- 
cess, where mass, momentum and energy are conserved, 
which leads to the build-up of correlations between the 
particles and gives rise to hydrodynamic interactions [115] . 

Three-dimensional periodic boundary conditions are 
considered for the simulation of shorter chains. Here, 
Lees-Edwards boundary conditions are applied to im- 
pose a shear flow [22]. A local Maxwellian thermostat 
is used to maintain the temperature of the fluid at the 
desired value [23]. A parallel MPC algorithm is exploited 
for systems of longer chains, which is based on a three- 
dimensional domain-decomposition approach Il6| . In such 
a system, shear flow is imposed by the opposite move- 
ment of two confining walls, and periodic boundary con- 
ditions are applied parallel to them. We impose no-slip 
boundary conditions at walls for both, fluid particles and 
monomers [T9l[25l[26] . 

Non-hydrodynamic simulations are performed by Brow- 
nian MPC, where each monomer independently performs 
stochastic collisions with a phantom particle which mim- 
ics a fluid element of size a 3 . In shear flow, the phantom- 
particle momentum is taken from a Maxwell-Boltzmann 
distribution with mean (p x ) = m (N c ) ^r y and variance 
(pp) = m(N c )k B T ((3 £ {x,y,z}), where (N c ) is the av- 
erage number of solvent particles per collision cell, r y the 
particle position along the gradient direction, and (p x ) the 
momentum along the flow direction [19l[27] . 

We employ the parameters a = 130°, h = O.lf, with 
f = sj ma 2 1 (kgT) (kg is Boltzmann's constant and T 




Fig. 1: Snapshot of a system of N p = 800 polymers of length 
N m — 250 for the Weissenberg number Wi c = 184 at the 
concentration c/c* = 2.77. For illustration, some chains are 
highlighted in red. Animations are provided as supplemen- 
tary material, low.mov shows the polymer dynamics under 
shear for the Weissenberg number Wi c = 18 and heigh. mov for 
Wi c = 184. 



is temperature), (N c ) = 10, M = m(N c ), I = a = a, 
kgT/e = 1, h/hp = 50, and the bond spring constant k = 
5 x 10 3 k B T/a 2 . The polymer lengths N m = 50 and 250 are 
considered in the concentration ranges c/c* = 0.16 — 2.08 
and 0.17 — 10.38, respectively. The corresponding over- 
lap concentrations are c* = 0.098/ 3 and c* = 0.029Z 3 , 
determined by their radii of gyration. In dilute solu- 
tion, the equilibrium end-to-end vector relaxation times 
are t = 6169f and 78330f [16]. The Brownian MPC sim- 
ulations for N m = 50 yield an approximately five times 
larger relaxation time than hydrodynamic MPC. 

Tumbling Dynamics. — The snapshot of a scmidi- 
lute solution, displayed in Fig. [1] indicates large confor- 
mational differences between the various polymers in flow. 
The average shape of an individual chain is illustrated 
in Fig. [2{a) by the monomer density distribution in the 
flow-gradient plane. Their conformational and rheological 
properties are discussed in detail in Ref. [16]. The large 
conformational variations are due to a continuous end- 
over-end tumbling motion [1 , 2^ , with stretched and coiled 
states as depicted in Fig. &b). 

The instantaneous shape of a polymer is characterized 
by the radius of gyration tensor Gppi ((3,(3' € {x,y,z}), 
which is defined as Gppi = ^iJHiAr^pAr^pi) /N m , where 
Ar^ is the position of monomer i in the center-of-mass 
reference frame of the polymer. To find a characteristic 
time for the tumbling dynamics, we determine the cross- 
correlation function 

c xy (t) = (G - (toK? ^ + t)) , (1) 

^(G' 2 x (t ))(G' 2 y (t )) 

for deviations from average stationary values G'pp(t) = 
Gpp(t) — (Gpp}. Figure [3] shows cross-correlation func- 
tions for several shear rates and concentrations. Each 
of the curves exhibits a deep minimum at time t + > 
and a pronounced peak at time t_ < 0, and decays to 
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Fig. 2: (a) Monomer density distribution in the flow-gradient 
plane for N p = 3000, N m = 250, i.e., c/c* = 10.38, and 
Wi c = 569. The contour lines for the densities 0.1 (outer) 
and 0.5 (inner) are highlighted to emphasize the non-ellipticity 
of the shape, (b) Illustration of polymer stretching (right) and 
recoiling (left). 9 is the angle between the end-to-end vector 
and its projection onto the flow-gradient plane and tp is the 
angle between this projection and the flow direction. 
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Fig. 3: Cross-correlation functions [Eq. {TJ] for a polymer of 
length N m — 250 and the concentrations c/c* = 2.77 (black 
solid line), c/c* = 10.38, (red dashed line), and c/c* = 0.35 
(blue solid line), corresponding to the Weissenberg numbers 
Wi c = 5520, Wi c = 5690, and Wi c = 2670, respectively, as 
well as N m = 50 with c/c* = 2.08 (green dashed-dotted line) 
and c/c* = 0.8 (green dashed line), the Weissenberg numbers 
are Wi c = 220 and Wi c = 233, respectively. Inset: Longest 
polymer relaxation times for N m = 50 (•) and N m = 250 (■). 



zero at large time-lags. Hence, the tumbling dynamics is 
not periodic, but cyclic. The latter has been questioned 
for tethered polymers [13j . The minimum at t + indicates 
that positive values of G^g are linked with negative ones 
of the orthogonal directions, i.e., polymer shrinkage in the 
y-direction is linked with its extension in x-dircction, and 
similarly, an extension in y-direction is linked to shrinkage 
in ^-direction. The maximum of C xy (t) reveals that posi- 
tive deviations G' xx are correlated with positive values G' yy 
at earlier times, or a collapsed state along the ^-direction 
{G' xx < 0) is correlated with a previous collapsed state in 
y-direction [2J. Hence, the time difference t + — £_ is re- 
lated to conformational changes that a polymer undergoes 
during tumbling. We therefore characterize tumbling by 
the time Tt = 2(t+ — t-). The factor two is introduced, 
because two non-equivalent conformations lead to a max- 
imum and a minimum, respectively, and will be (more or 
less) assumed during a cycle. 

As shown in Fig. [3] the positions of the peaks and min- 
ima are rather close for equal Weissenberg numbers, when 
the lag-time is scaled by the relaxation time t(c). Hence, 
the tumbling times exhibit a strong concentration depen- 
dence due to the strong concentration dependence of the 
relaxation times, which is shown in the inset of Fig. 3 [16| . 

Normalized tumbling frequencies / = r(c)/r t , with 
tumbling times extracted from the correlation functions 
and scaled by the corresponding relaxation times r(c), are 
presented in Fig. [4] for a wide range of shear rates and 
concentrations. For comparison, the theoretical predic- 
tion for a polymer in dilute solution is presented as well 
[12j[2S] . The results are in excellent agreement, when the 
Weissenberg number Wi* of the theoretical model is iden- 



tified with Wi* = Wi c /2, where the factor two accounts 
for the approximately twice larger relaxation time of the 
theoretical model. The short chain results clearly show 
the crossover from unity, assumed in the limit j — > 0, to 
the asymptotic dependence ~ Wi^ 3 at high shear rates. 
We obtain a chain-length dependence in close agreement 
with the theoretical prediction. More importantly, we find 
a slight and gradual shift of / to larger values with in- 
creasing concentration at a given Wi c , until a saturation is 
reached in the semidilute regime c/c* > 1. This is seen for 
the two largest concentrations for N m = 50 and the three 
largest ones for N m = 250. As a consequence, the poly- 
mers exhibit a universal behavior, both in dilute (c <C c*) 
as well as in semidilute solution as function of Wi c , with 
the same power-law dependence on Wi c (for Wi c > 1). 

Theory [12], simulations [10], and experiments [5] sug- 
gest that tumbling is an aperiodic process with an expo- 
nential distribution of intervals between tumbling events, 
i.e., Pt(t) ~ cxp(— i/r t e ), where r t e is defined as the char- 
acteristic tumbling time. By calculating the distribution 
functions of times between successive gradient-vorticity 
plane crossings of the end-to-end vector, we determined 
the tumbling times r t e presented in the inset of FiglH Their 
dependence on Weissenberg number and concentration is 
in perfect accord with that obtained for r t ; the absolute 
values are only approximately 20% smaller. 

Alternatively, relaxation times under shear flow can be 
obtained by the end-to-end vector auto-correlation func- 
tion (R (t)R (O)), where Rp = r N ,p - r hf3 [H1H21E9]. 
Similar to the results presented in Ref. [29], we find 
a damped oscillatory behavior for Wi c > 1. A fit of 
git) = e~*/ r ''[cos(wi) + asin(wi)] to the correlation func- 
tions along the flow and gradient directions, yields, within 
the accuracy of the simulations, relaxation times r r (7) and 
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Fig. 4: Normalized tumbling frequencies / = r(c)/r t obtained 
from cross-correlation functions. Open symbols correspond to 
the polymer length N m = 50 for c/c* = 0.16 (o), c/c* = 0.4 
(V), c/c* = 0.8 (A), c/c* = 1.6 (0), and c/c* = 2.08 (□). 
Filled symbols indicate results for N m — 250 and c/c* = 0.35 
(•), c/c* = 0.69 (►), c/c* = 1.38 (▼), c/c* = 2.77 (♦), c/c* = 
5.19 (+), c/c* = 10.38 (■). The lines present the theoretical 
predictions [12] [28]. Inset: Relaxation times r t e for N m = 50 
and various concentrations. 



normalized tumbling frequencies / = T(c)/r r equal to the 
values presented in Fig. 2J The parameter u of g(t) is in- 
dependent of shear rate for the considered systems. This 
is in contrast to results of Ref. [55] , where very short poly- 
mers (N m = 10) have been considered only. 

The simulation results for the tumbling time lead to the 
following conclusions, (i) The correlation function ([1} can 
be used to obtain a characteristic time Tt for the tum- 
bling motion (see also Ref. [3]). (ii) The non-equilibrium 
end-to-end distance relaxation time r r along the flow di- 
rection is equal to the tumbling time r t extracted from 
the correlation function (JT|). Hence, the tumbling time r t 
is equal to the non-equilibrium relaxation time r r , as also 
predicted in Refs. [H1Q3]. (iii) There is a weak chain- 
length dependence of the tumbling time, (iv) Screening of 
hydrodynamic interactions in semidilute solutions leads to 
a (small) increase of the normalized tumbling frequencies. 
The screening aspects are discussed in more detail in then 
next section. 

Angular Probability Distribution Functions. 

Further insight into the tumbling and orientational behav- 
ior of polymers is gained by the probability distribution 
functions (PDFs) P((p) and P{6) for the orientation angles 
<p and 6 (cf. Fig. (Kb)) [gmOSQl]. For a dilute solution, 
P((p) is shown in Fig. [SJ together with theoretical lines 
obtained from Ref. [12] (note, Wi* = Wi c /2). Evidently, 
the simulation results agreement well with the analytical 
approach, as is expected for a dilute solution in which the 
intermolecular interactions are irrelevant. P(ip) exhibits a 
significant shear-rate dependence. At zero shear, no angle 
is preferred. An increasing shear rate leads to the ap- 
pearance of a peak, which shifts to smaller values with 
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Fig. 5: Probability distribution functions of the angle tp for a 
dilute solution with c/c* = 0.16 of polymers of length N m = 50 
for Wi c = 616.9 (■), Wi c = 61.7 (o), and Wi c = 12.3 (A). 
Inset: Distribution functions for N m — 250, c/c* = 0.17, and 
Wi c = 2350 (o), Wi c = 235 (x), Wi c = 23.5 (□). The solid 
lines are theoretical results for Weissenberg numbers Wi* = 
Wi c /2. 



increasing j and, at the same time, the width Atp of P((f) 
decreases. 

In a semidilute solution, the universality observed for 
the tumbling time is also reflected in P((f). Figure [5] dis- 
plays distribution functions for various concentrations and 
Weissenberg numbers Wi c . For every Weissenberg number 
distributions are compared for three concentrations. Evi- 
dently, the distributions are independent of concentration 
for the considered Weissenberg numbers. 

However, we observe a clear concentration dependence, 
when we compare distributions from dilute and semidilute 
solutions. The inset of Fig. [6] displays distributions for 
the concentrations c/c* — 0.35, 5.19 and the Weissenberg 
numbers Wi c « 267, 2670 and 270, 2700, respectively. 
Clearly, the increase in concentration from a dilute solu- 
tion beyond the overlap concentration leads to a broaden- 
ing of the distribution function. 

Interestingly, the value tp m at the peak of the distribu- 
tion function is independent of concentration at a given 
Wi c . Figure [7] displays tan(2</j m ) as function of Wi c for 
various concentrations of the two studied chain lengths. 
The simulation data are consistent with the analytical pre- 
dictions for both chain lengths [12l[28]. The results show 
that a universal behavior is obtained for the various con- 
centration at a given N m . In the asymptotic limit of high 
shear rates, the dependence 

tan(2v? m ) ~ (Wie)" 1 / 3 (2) 

is obtained, which has also been found in Ref. [16] for the 
alignment angle determined from the gyration tensor. In 
the limit of Wi c — > 0, theory predicts tan(2y> m ) ~ Wi^ 1 , 
whereas the simulations yield tan(2cp m ) ~ Wi~ ' 8 for the 
considered range of Wi c , which might be due to excluded 
volume interactions not taken into account in the theo- 
retical calculations. The inset of Fig. [7] displays A(f, the 



p-4 



Tumbling Dynamics of Polymer in Semidilutc Solution 




Fig. 6: Probability distribution functions P(ip) of polymers 
of length N m = 250 for c/c* = 2.77 with Wi c = 5520 (■), 
Wi c = 552 (♦), and Wi c = 55.2 (a), for c/c* = 5.19 with 
Wi c = 5423 (o), Wi c = 542.3 (□), and Wi c = 54.23 (+), as 
well as c/c* = 10.38 with Wi c = 5691 (solid line), Wi c = 
569.1 (dashed line), Wi c = 56.91 (dashed-dotted line). Inset: 
Distribution functions for c/c* = 0.35 and Wi c = 2670 (black 
solid line), Wi c = 267 (red solid line) as well as for c/c* = 5.19 
with Wi c = 2700 (black dashed line) and Wi c = 270 (red 
dashed line). 



full width at half maximum. For dilute solutions, Aip 
agrees with the prediction of the theoretical model, which 
yields the asymptotic dependence Aip ~ Wi^ 1 / 3 in the 
limit Wi c — > oo. The widths of the distributions are larger 
for semidilute solutions, but the asymptotic Weissenberg 
number dependence seems to be the same. At c/c* > 1 a 
universal curve is adopted, as already pointed out above 
in connection with Fig [51 

Probability distribution functions of the angle 8 are dis- 
played in Fig. [5] for various Weissenberg numbers. Theo- 
retical calculations for dilute solutions predict a crossover 
from a Gaussian shape of the distribution function to a 
power-law decay according to P{6) ~ 9~ 2 with increas- 
ing shear rate, within a certain range of angles, which is 
confirmed by the simulations. The dependence of P(6) 
on concentration is in accord with that of P(<p). In the 
semidilute regime, P{9) is independent of concentration 
for a given Wi c , while a comparison of distributions of 
semidilute and dilute solutions yields a broadening for the 
semidilute case. 

We attribute the concentration independence of the 
tumbling time and the probability distribution functions 
for c/c* > 1 to screening of hydrodynamic interactions. 
To confirm our hypothesis, we performed Brownian MPC 
simulations for dilute and semidilutc solutions. Using sim- 
ilar Weissenberg numbers, we find, within the accuracy of 
the simulations, identical distribution functions P(<p) for 
both cases. Moreover, the distributions agree with those 
of semidilute systems of the same concentration and Weis- 
senberg number in the presence of hydrodynamic interac- 
tions. Hence, the differences between distribution func- 
tions at low and high concentrations, as displayed in the 
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Fig. 7: Maximum angle ip m as a function of Wi c for N m = 250 
and c/c* = 0.17 (o), 0.35 (▼), 2.77 (♦), 5.19 (a), 10.38 (■), as 
well as Nm = 50 for c/c* = 0.16 (►), 1.6 (x), and 2.08 (□). 
Inset: Widths Aip of the distribution functions. For both, the 
lines are theoretical predictions. 



insets of Figs. 6 and 8, arc due to hydrodynamic interac- 
tions. In dilute solutions, hydrodynamic interactions are 
present, whereas in systems with c > c* , hydrodynamic 
interactions are screened. Naturally, at larger concentra- 
tions friction is higher. This aspect is captured in the 
relaxation time t(c), which increases considerably with 
concentration [16J . 

The fact that the tumbling frequencies and the widths 
of the distribution functions are larger in semidilute solu- 
tions, i.e., when hydrodynamic interactions are screened, 
might be explained as follows, (i) The observed broad- 
ening of the distribution function P{ip) in semidilute so- 
lution implies that hydrodynamic interactions favor poly- 
mer alignment and lead to a faster dynamics during the 
collapse and stretching part of the tumbling motion, (ii) 
At the same Weissenberg number, the shear rate of a non- 
draining polymer is larger than that of a free-draining one, 
due to differences in equilibrium relaxation times, i.e., for a 
coiled conformation. As a consequence, the effective Weis- 
senberg number Wir = 7T#, where tr is the rotational 
relaxation time in the stretched rodlike conformation, of 
the non-draining polymer is larger than that of the free- 
draining one. This could explain the larger probability of 
angles in the vicinity of ip m for non-draining polymers as 
well as their faster collapse dynamics. Overall, the tum- 
bling time is larger in a non-draining system. 

The broadening of the distribution functions with in- 
creasing concentration or screening of hydrodynamic inter- 
actions is not captured by standard theories employing the 
preaveraging approximation [j"2 l 1171128] . Here, hydrody- 
namic interactions arc included in the relaxation times and 
hence the Weissenberg number only. Additional, "higher 
order effects" arc neglected. Therefore, one might expect 
that the theoretical description would reproduce results 
of simulations without hydrodynamic interactions, in con- 
trast, the model calculations rather reproduce the Simula- 
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Fig. 8: Probability distribution functions P(ff) of polymers of 
length N m = 250 for c/c* = 2.77 with Wi c = 5520 (■), 552 
(♦), 55.2 (a), and 0.552 (o), as well as c/c* = 10.38 with 
Wi c = 5691 (solid line), 569 (dashed line), 57 (dashed-dotted 
line), and 0.57 (thin line). The slope of the short solid lines 
is —2. Inset: P(0)s of a dilute solution with c/c* = 0.35 for 
Wi c = 2670 (black solid line) and Wi c = 267 (red solid line) 
and a semidilute solution with c/c* = 5.19 for Wi c = 2700 
(black dashed line) and Wi c = 270 (red dashed line) . 



tion data for systems with hydro-dynamic interactions. 

V. Conclusions. — Wc have found in Ref. [TB] that 
in shear flow the stationary-state conformational poly- 
mer properties arc independent of concentration when ex- 
pressed in terms of the Weissenbcrg number Wi c = 7t(c). 
This is remarkable, since the longest polymer relaxation 
time t(c) increases significantly with concentration and 
indicates that an effective local friction determines the 
stationary-state properties. 

Here, we have analyzed dynamical properties — 
orientational distribution functions and tumbling times — 
of semidilute polymer solutions in shear flow and have 
found that they depend on concentration (in excess of 
t(c)), a dependence which we attribute to screening of 
hydrodynamic interactions in semidilute solution. Com- 
pared to the dilute case, such a screening causes a broad- 
ening of orientational angle distribution functions and an 
increasing ratio / = T t /r(c) at the same Weissenbcrg 
number Wi c in semidilute solution. The effect itself is 
small (/(c = 0)//(c > c*) « 1.3 at Wi c = 10 3 ). More 
importantly, the same asymptotic dependencies are ob- 
tained as function of the Weissenberg number Wi c in di- 
lute and semidilute solutions. This explains the previously 
obtained agreement of power spectral densities obtained 
from free-draining and non-draining computer simulations 
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